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We investigate the entanglement entropy of a block of L sites in quasifree translation-invariant 
spin chains concentrating on the effect of reflection symmetry breaking. The majorana two-point 
functions corresponding to the Jordan- Wigner transformed fermionic modes are determined in the 
most general case; from these it follows that reflection symmetry in the ground state can only be 
broken if the model is quantum critical. The large L asymptotics of the entropy is calculated ana- 
lytically for general gauge-invariant models, which has, until now, been done only for the reflection 
symmetric sector. Analytical results are also derived for certain non-gauge- invariant models, e.g., 
for the Ising model with Dzyaloshinskii-Moriya interaction. We also study numerically finite chains 
of length with a non-reflection-symmetric Hamiltonian and report that the reflection symmetry of 
the entropy of the first L spins is violated but the refiection-symmetric Calabrese-Cardy formula is 
recovered asymptotically. Furthermore, for non-critical reflection-symmetry-breaking Hamiltonians 
we find an anomaly in the behavior of the "saturation entropy" as we approach the critical line. The 
paper also provides a concise but extensive review of the block entropy asymptotics in translation 
invariant quasifree spin chains with an analysis of the nearest neighbor case and the enumeration of 
the yet unsolved parts of the quasifree landscape. 



I. INTRODUCTION 

Understanding the entanglement properties of systems 
with many degrees of freedom, such as quantum spin 
chains, has been one of the main recent research top- 
ics connecting quantum information theory and con- 
densed matter physics [1-5]. Huge amount of results 
have been accumulated about translation-invariant sys- 
tems. However, the results almost exclusively correspond 
to reflection symmetric systems, despite the fact that 
models violating reflection invariance play a prominent 
role in many-body theory, e.g., in describing interactions 
of Dzyaloshinskii-Moriya type or non-equilibrium steady 
states. 

Considering a subsystem 5 of a system, which is in 
a pure state, the entanglement between the subsystem 
and its environment is characterised by the von Neumann 
entropy 

S{ps) '■= -Tr(p5 lnp5) , 

where ps denotes the density matrix of the subsystem. 
In the case of infinite one-dimensional critical chains, this 
entanglement entropy belonging to a block of L contigu- 
ous spins was shown to grow asymptotically as [1, 2] 

SL = ^\nL + k, (1) 

where c is the conformal charge of its universality class 
and A; is a non-universal constant. For non-critical chains 
the asymptotics of the entanglement entropy is bounded. 
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This saturation value of the entropy diverges as one ap- 
proaches the critical point: it increases as [2] 

^,,t = ^lne + A:' , (2) 

where ^ is the correlation length. In the case of finite 
chains (with open boundary conditions) consisting of N 
spins, the conformal field theoretic prediction for the en- 
tanglement entropy of the first L spins (at criticality) is 
[2, 6, 7] 

, c /2A^ ttL\ k 
5(L,iV) = -ln(^— sin— j -l-lng-f- , (3) 

where In g is the boundary entropy introduced by Affleck 
and Ludwig [8]. 

In this paper, we will study the asymptotics of the en- 
tanglement entropy in chains with broken rcfiection sym- 
metry. We consider quasifree models (with finite range 
coupling) : their Hamiltonian can be mapped to quadratic 
fermionic chains by the Jordan- Wigner transformation 
[32] 




Throughout the paper we will assume either open bound- 
ary conditions or "fermionic" periodic boundary condi- 
tions {bi — bi+N) [33]. The requirement of translation- 
invariance implies that A and B are Toeplitz matri- 
ces {Ai+n,j+n = and B^+^.j+n = Bij for any 
n € N), hermiticity of H implies that ^ is a (possibly 
complex) hermitian matrix, and B is (a possibly com- 
plex) anti-symmetric matrix. Finite-ranged interaction 
means that there exists a positive integer no such that 
Ao,i = Bqj = if / > ng. Such a spin-chain Hamiltonian 
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is not invariant with respect to the reflection transforma- 
tion R{iy'^) = cr° j (a = X, y, z), iff A is not a real matrix. 
(One might think tliat the term {hibj — hjbi), with i > j 
also breaks the translation invariance of the spin chain, 
but a short calculation shows that its image under the 
Jordan- Wigner transformation is the following reflection- 
invariant term (Ji^'~2i+^<Jl(Jj + <^J^n=j-i<<^i)- One 
of the most studied quantum spin chain with broken re- 
flection symmetry is the Ising model with transverse mag- 
netic field and Dzyaloshinskii-Moriya (DM) interaction 
(in the 2;-direction) [10-13]: 

N 

H = Y1 «+i + f^'^i + D{crfa^+, - a^af^,) , (5) 

i=l 

Another type of model that has been studied extensively 
in the literature is the model 

H = Y.[j{afat^,+ayaf^,)+ (6) 

i=l 

whose ground states are used to describe the energy cur- 
rent carrying eigenstates of the XX model [14-16]. Cer- 
tain non-reflection invariant quasifree states also appear 
as invariant states of reflection-invariant quantum cellu- 
lar automata [17]. 

The entanglement entropy asymptotics of the models 
given by Eq. (4) has been studied by many authors [18- 
22]. The main analytic tool for tackling this problem was 
expressing the entropy in terms of the determinant of a 
Toeplitz matrix, applied first by Jin and Korepin [18]. 
Until now the most general results have been achieved 
by Keating and Mczzadri [19], who gave a general ana- 
lytic expression for the entropy asymptotics when A is 
real and B = 0, and by Its, Mo, and Mezzadri [22], who 
gave an analytic (although less explicit) expression even 
for the case of general (finite-ranged) real A and B ma- 
trices, while certain results about the d dimensional case 
can be found in [23]. However, none of these studies con- 
cerned reflection symmetry breaking cases, i.e., when A 
is complex. 

We will generalize the above mention results by deriv- 
ing an analytic expression for the general gauge-invariant 
case (i.e., when ^ is a "general" complex Hcrmitian 
finite-ranged Toeplitz matrix, while B = 0). This in- 
cludes, as a particular case, the model described in 
Eq. (6). Moreover, we will also introduce a multitude 
of transformations between models of Eq. (4) , which al- 
lows for deriving analytic expressions for cases with non- 
vanishing B. A remarkable result that we obtained is 
that for these " quasifree" models reflection invariance can 
only be broken in the ground state if the model is criti- 
cal. If the model is non-critical the ground state of the 
model does not change if we replace Aij with I{je{Aij) in 
the Hamiltonian. From this, as we will show, it follows 
that scaling in Eq. (2) may be violated. However, we will 



discuss how we can reinterpret this equation to keep its 
validity. Furthermore, we will present numerical results 
in non-reflcction-symmetric spin chains providing an ex- 
ample of broken reflection symmetry in the finite size 
scaling of the entropy S{L, N) ^ S{N — L, N) breaking 
the symmetry of Eq. (3), but we will see that that this 
deviation goes to zero as we increase the system size. 

The paper is structured as follows. In Section II we cal- 
culate the majorana two-point functions of these general 
(finite-ranged) quasifree models and recapitulate how one 
can obtain the entanglement entropy from the two-point 
functions. The results already known about the entan- 
glement asymptotics of certain types of quasifree models 
arc collected in Section III. We derive an analytic formula 
for the entanglement entropy for general gauge-invariant 
models in Section IV, whereas in Section V we show how 
we can extend our results for certain types of non-gauge- 
invariant models too. Section VI is an application of the 
above to models with nearest neighbor interactions, while 
in Section VII we discuss how some of our analytic and 
numerical results conflict with the formulas (2) and (3) 
and how we can "resolve" this discrepancy. Finally, Sec- 
tion VIII is devoted to the summary and the remaining 
open questions. 

II. TWO-POINT FUNCTION OF THE 
MAJORANA OPERATORS AND 
ENTANGLEMENT ENTROPY 

The entanglement entropy asymptotics of the models 
described by the quadratic Hamiltonians in Eq. (4) can 
be calculated from the ground-state expectation values 
(mfcm/), where m„'s denote the so-called majorana op- 
erators defined as 

m2n = i{bn-b'l), m2„-i = bn + bl . (7) 

In this section we will first derive these majorana two- 
point functions in terms of the matrices A and B that 
define the Hamiltonian Eq. (4). Thcin wc; describe how 
to calculate (in this quasifree setting) the entanglement 
entropy alone from two-point functions, and finally we 
recapitulate the "determinant trick" of Jin and Korepin, 
which will allow us later to obtain analytical results. 

A. The majorana two-point functions 

Let us fix our conventions used in the calculation. We 
will consider the "fermionic" periodic boundary condi- 
tion: bi = bi+N- The Fourier and inverse transforms of 
the one-particle annihilation operators read 

^'^ = ^l.^^P[-^)bn (8) 

T 1 f 2TTink\ ~ 
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the summation runs in the set of integers [■^^j^j'^^j^ 

( [~ T ' T ~ -*■] ) ^'^^ ^ '-"^'^ (even) and the transform of 
the one-particle creation operators are to be computed 

by means of taking the adjoint of the above formulae. 
For Toeplitz matrices we define the Fourier transform as 



] where = {Ak + A^k)'^ + ABkBl, and the two remain- 



Xk = ^ exp r - ^^^^ jXo,n 



0,n 



u \ 



N J 

( 2Trink\ 

IT) 



Xk 



(10) 

(11) 



here Xo,„ stands for either Ao,„ or Bo.n, and the sum- 
mation again runs in the set of integers [-^^^,"^^2^] 
([— ^, ^ — l] ) for N odd (even). Using these definitions, 
the Hamiltonian (4) can be written as 



H 



E 

k 



(12) 



To bring this Hamiltonian into a diagonal form H = 
Aft c^Cfc, {Ak G M), one performs a Bogoliubov trans- 
formation 



Ck = Oikh + Pkb^-k ak,Pk&'C 
where the coefficients ak,Pk have to satisfy 

akP-k + PkOt-k = 
|afc|' + |/3fc|' = 1 , 



(13) 

(14) 
(15) 

so that the canonical anticommutation relations 
{cfe,c^,} = dkk' are satisfied. The consistency conditions 
for the commutator [cfe ,H]=Ak Ck give 



-Ak Bl 
Bk A-k 



Oik 

13k 



Ak 



ak 



One readily extracts the one-particle spectrum 

A_k -Ak + ^{Ak + A_k? + ^BkBl 
At = 



(16) 



(17) 



having taken the relations A*^. = Ak, B_k = —Bk (which 
are direct consequence of Aj ^ = Aij,Bj^i = —Bij) into 
account. The ground state correlations for the two-point 
functions of the new Fermi operators read 

and all other correlations vanish. Now, using the inverse 
of (13), bk — c(*kCk + f3-k c^-fe, one can compute the cor- 
relations among the Fourier components bk, bj., and sub- 
stituting the solution of (16) for ak, Pk we arrive at 

1 



- 2Trik{j - I) Bk 

k 



Afc ^ A-fc 
|Afe| |A_fe| 



V- 2mk(j - 1) 

> exp X 

N ^ N 

k 



{b]bi) 4 



(18) 



-fc \ I A-k+Ak ( Afc I A-fc ^ 



ing two-point functions and {bjbj) can be calcu- 

lated directly from the above equations. Ultimately, we 
would like to have a linear combination of the above, the 
two point functions of the self-adjoint majorana opera- 
tors defined in Eq. (7). Before writing down the final 
result, let us introduce some notations. We will use the 
combinations 



Al = A-k + Ak, Al = A-k - Ak 
Bk = Bk + Bl , i B^ ^ Bk - Bl 



and the step functions 
A-fc 



PkA 



|A-fc| 

Note, that Al,Al,Bl, B^ G 



Afc A-fc 
lAfcl + |A_fc| 

and Afc 



(19) 
(20) 



(21) 



(B|)^ -I- {B^y. We now take the thermodynamic limit 
{N oo) and write the final result in a manner usually 
adopted in the literature 



{m^m,i) = Sji + iCji , 
where the matrix C has the following structure 



(22) 



Ho 




n-2 


n-3 


Hi 


no 


n_i 


n_2 


n2 


Hi 


Ho 


n_i 


ns 


n2 


Hi 


Ho 



c = 



The 11; 's are 2x2 block entries that read 



(23) 



-iie . 



ui^^ Jdee- 



^/m 



(24) 



where 9 — 2'Kik/N , so all Fourier series become functions 
on the circle [0, 2tt] in the limit. This type of matrix C 
is called block- Toeplitz and the matrix argument in (24) 
of the integral ip : ^ M2(C) is called its symbol. 

The n-point majorana function can be obtained from 
the two-point functions by the Wick rule [9]: 

fc 

Esgn(7r) n ("''7r(2i-i)m^(2;)) . 

where the sum runs over all pairings of {1,2,..., 2k}, i.e., 
over all permutations of the 2k elements which satisfy 
tt{21 - 1) < n{2l) for / € {1, 2, . . . , fc} and 7r(2Z - 1) < 
tt{21 + 1) for ; G {1,2,...A;- 1}. 
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Before coming to the calculation of the entropy, let 
us analyse the obtained result. The one particle spec- 
trum (17) has the form of a sum of a reflection invari- 
ant y/A{6)/2 and a non-invariant term A'^{6)/2 (Note 
that the real space reflection n — > —n corresponds to 
the Fourier space one k — > — fc as follows from the 
Fourier transform (8)). The symbol (24) characteriz- 
ing the correlation matrix {miinj) has a dependence on 
the non-reflection invariant part of the spectrum only 
via Me = (A(6')/|A(6i)| - A(-6')/|A(-^)|)/2. This term, 
however, vanishes identically unless A(^?o) = at some 
6*0 G [0,27r]. In other words, non-critical quasifree sys- 
tems never break reflection invariance [34]. We will dis- 
cuss some implications of this important fact in Sections 
VI and VII. 



-1 



c 

■r\r(e) 



r=e/2 



FIG. 1: The integration contour encloses a domain, where 
e(l -|- e. A) defined in (26) is analytic and contains all singular 
values of the matrix Ci, . 



this form.) The entropy can now be easily calculated. It 
can be written in terms of the function 



e{x, v) 



X + v ^ x + u X — u ^ X — u 
m in 



(26) 



B. Calculation of the entanglement entropy from 
the two-point functions 

Restricting the ground state to a subsystem consisting 
of L consecutive sites one obtains a mixed state. Let us 
restrict the matrix C defined in Eq. (22) (which describes 
the two-point majorana correlations) to L consecutive 
modes, that is to a 2X x 2L submatrix 



as 



Sl = S{pl) = ^e(l,z/j) . 



(27) 



The trick [18] to obtain the asymptotics of the entan- 
glement as the size of the block grows is computing the 
determinant 



Cl = 



( Ho n_i 

Hi Ho 



n_L+i\ 
n_L+2 

Ho } 



(25) 



where the 11; 's are 2x2 matrices given by Eq. (24). Let 

us denote by W the orthogonal matrix, the adjoint action 
of which brings the antisymmetric real matrix Cl into its 
canonical form, i.e., for (fl'i,)ij = Z^fe^Io Wik{CL)kiWji, 
we have 



Hr 



fe=l 



1 
-1 



where Vk G [0, 1] (fc = 1, 2, . . . , L) are the singular values 
of Cl- (Due to the fact that Cl is antisymmetric, the 
degeneracy of its singular values is always an even num- 
ber, that is why we label them only from 1 to L.) The 
density matrix corresponding to the restricted state can 
be written as [35] 



PL 



n 



l + Vj 1 + im2j-iiri2j _|_ 1 



Vj 1 — im2j-im2j 



where fhj 



J2lto^ Wjkmk for all j = 0, 1, ... 2L - 1. 
(Actually, translation invariance is not used here, the 

density matrix of any quasifree state, i.e., of any state 
for which the Wick expansion applies, can be written in 



L-l 

Dl{X) = det {iXI + Cl) = (-1)^ J] ' ' (^8) 

3=0 

and exploiting the residue theorem by writing down the 
following integral 



lim 



1 



where the contour T{e) is shown in Fig. 1. That contour 
encircles all eigenvalues of C^, but bounds a region, in 
which e(l -|- e. A) is analytic. Hence, the main task in all 
cases is to compute the determinant of the block- Toeplitz 
matrix matrix 



(30) 



Finally, we should mention that in the gauge-invariant 
case (i.e., when B — Q) there is an easier method for the 
calculation of the entanglement entropy. In this case, 
as can be seen from Eq. (18), the non-gauge-invariant 
two-point functions vanish {{bjbk) = = 0). If 

we restrict the state to L consecutive sites, and de- 
note by Ml the corresponding restriction of the ma- 
trix Mij = {blbj) and by Ul the (not necessarily real) 
unitary, the adjoint action of which diagonalizes Ml 

iY.j,k=i Uij{ML)jk Uik = XiSii), then the density matrix 
of tne restricted state reads 



PL = l[{Xi4c, + {l-Xi)c,4) 
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where Ci = J2j Uijbj and Xi € [0, 1]. Hence the entropy 
of the restricted state is given by 

L 

Sl = -TrpLlnpL = -J2 i^i In Ai + (1 - Xi) ln(l - A^)) 

(31) 



III. SUMMARY OF PREVIOUSLY KNOWN 
CASES 

In this section we shortly recapitulate what has been 
previously known about the entanglement entropy for 
quasifree models. 



A. Gauge and reflection invariance 

In the case of gauge- and reflection-invariant quasifree 
models, the matrix B is zero while A is real, which im- 
plies B'ie) = B^ie) = A" (9) = and the symbol of the 
majorana two-point functions (24) reduces to 





Hence, Cl can be factorized as 



Cl{X) 





-1 



\A-{8)\ 





' Gl(A) 



(32) 



where Gz, ( A) is the restriction of the Toeplitz matrix with 
scalar symbol g{e) = A''{e)/\A%9)\ to an L x X block on 
the diagonal. From this it follows that det(CL) = = 
(-1)^ det(A/ + GL) det(A/-Gi,). To extract the entropy 
asymptotics, one only needs to calculate the (L — >■ oo 
asymptotics of the) determinant of (A/ ± Gl) using the 
Fisher-Hartwig theorem, and then use the residue theo- 
rem as described in Section II B. This was done by Keat- 
ing and Mezzadri [19, 20]: they obtained the following 
result: Let there be number R/2 zeros of A" (6*) denoted 
by {r = l, - ■ ■ , R/2) in the semi-circle [0, n] (implying 
another R/2 zeros in the other semi-circle [— 7r,0]: —0^, 
r — 1, . . . , R/2). Then the entanglement entropy asymp- 
totics is given by 



where 



54p^) = I In L I (In 2) /3 , 



R/2 



(33) 



K 



- Yl (-l)''+^ln 



l<s<r<R/2 



1 _ gt{0r~e,) 



1 _ e^ie^+e.) 



where je = ■■■ is Euler's constant, and I3 = 0.0221603..., 
independent oi A{6) (consult [18] for its derivation). 



B. Reflection invariance and real Bij 

The other case that has already been discussed in the 

literature is the case when both matrix A and B are real, 
i.e., when A°-{9) = B'^{9) = 0. In this case, the symbol 
reads 





\A^{e)+iB<'le)\ 



A'{e)-iB'(e) 




(34) 



Here the idea, invented for the XY model in [21] and 

generalized for the present case in [22], is to extend 
the domain of : S*^ — )■ M2(C) to the complex plane 
and use a theorem of Widom [26], which yields a for- 
mula of the block Toeplitz determinant at hand ex- 
pressed in terms of Wiener-Hopf factors of the symbol: 
ip{z) — U+{z)U-{z) = V-{z)V+{z), where the matrices 
U+,V+ {U-,V-) are analytic inside (outside) the unit 
circle. The factorization resides on the fact, that due to 
the assumption of finite range interaction, the functions 
A{z) = A{exp{-i9)) = A{9), and B{z) = B{exp{-i9)) = 
B{9) are Laurent polynomials. One writes 



A'^iz) - B'^iz) _ q{z) 



q{z) 



\A{z)-B{z)\- \q{z)\ \lq{l/z) \ 



2no 



with Zi being the roots of the polynomial p(2;) = 2"°(j(z), 
where no is the range of the coupling (defined after (4)). 
Note, that the equality in the middle is a choice of an- 
alytic continuation as q{z)* = q{l/z) holds on the unit 
circle (as is obvious from the general form g(exp —iO) = 
A* (9) — (9) . The non-analytic behaviour of the above 
rational function is then the only thing that has to be 
taken care of and the factorization is done with the help 
of theta functions living on the hyperelliptic surface of 
genus no given by 



w 



2no 

n(- 



Zi){l 



ZZi 



(35) 



The XF-model has no = 1 thus the underlying Riemann 
surface is a torus, while for for general finite ranged cou- 
plings q{z) can be any degree no Laurent polynomial, 
which satisfies q{z)* = q{l/z) on the unit circle. The 
result (Theorem 3. of [22]) for the logarithmic derivative 
reads 



dX 



2\L I 
1^ + 



dU+(z)jj-l 



dz 



U-'iz) + V-\z 



dz 



)g-\z)_ 



dz 
(36) 

and the difference (rhs.— lhs.)< Cp ^ where the constant 
p satisfies 1 < p < min{|Ai| : |Ai| > 1} (the complex 
numbers A, are the roots of p{z) and their reciprocals). 
The saturation entropy is given by 



S{pa) = 2 



' 9(^(A)^ + §)9(/3(A)^-i) 
©Mi) 
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This formula depends on the surface (35) via the theta 
functions (which are uniquely defined by some quasi- 
periodicity properties along non-contractible curves on 
the surface); their definition and that of their arguments 
will be omitted here (see [22]). We only remark that it 
is exactly at criticality, when the above surface becomes 
degenerate and the formula diverges. 



IV. GAUGE INVARIANT MODELS IN 
GENERAL 

The reason why one could give an explicit formula for 
the entropy asymptotics in the reflection and gauge in- 
variant case (when A" (6*) = B^iO) = B'{e) = 0) and a 
less explicit one in the case when A'^{9) = B^ii)) = was 
that the structure of the symbol (fi{6) was considerably 
simplified in both cases. 

In the general quasifree case it is hard to find the 
Wiener-Hopf factorization of the symbol, since there is 
no identically zero entry of in the matrix function (p{9). 
This is true even in the restricted case of gauge invari- 
ant (but not reflection invariant) models. However, as 
we will show in this section, one can circumvent this 
problem in this restricted case. We have seen in sec- 
tion II B that we can extract the entropy also from the 
correlation matrix = (&i&])|ij=i..L: it is given by 

Sl = - J2i=i i^i 111 + (1 - ^i) ln(l - Aj)), where Aj 
are the eigenvalues of the matrix C^. Now, we can use 
the contour integral trick again with a small alteration 
and write the entropy as 



5i, = lim-^-/ e(l-|-e,A) 



rfln£)^(A) 
dX 



(37) 



where D^(A) = dctC£(A) = det(A/ - (2C^ - /)), the 
function e{x, A) and the contour F were defined in Section 
II B. Hence the situation is analogous to section III A 
except that XI — Gl is replaced by XI — (2C^ — J), which 
is also a Toeplitz matrix, but its symbol 



X + l-2c'{e) = X- 



m_ 



(38) 



is not necessarily symmetric {Aij ^ Aj^i implies c'{9) ^ 
c.{—6)). Now we can use the Fisher- Hartwig conjecture 
[25]: Suppose that the symbol p(^) : 5^ C of a Toeplitz 
matrix has the following form 



P{O) = m]l*Pr,eA0)uc.,eA^) (39) 



r=l 



with 



Ua^,eA&) = (2-2cos{0-9r)y\ Rear > 



where the function xp : ^ C is smooth, non-vanishing 
and has zero winding number. Then the L — >■ oo asymp- 
totic formula for the determinant reads 

detPi = (^exp^^y^ liii;{e)d9^^ x 

||^^L«^^?jf[V',{aO,{A},{ei}], 

where 

£[ij,{ai},{(3i},{9i}] = £[ij] X 

R 

X JJ Gb(1 + ai + A)Gs(l + ai- Pi)/GB{1 + 2ai) 



i=l 
R 



-ai-f)i 



-ai+i3i 



xn(V-(exp(i^i))) (v+(exp(-i^i)) 

x[] h - exp(i(^i - ^, )j 

and the so-called Barnes function is defined by 
Gb{1 + z) = (2^)-/2e-(^+i)^/2-^-^'/2 X 

■DC' 

[|{(l + 2/n)"e-"+^'/(2n)|_ 

n=l 

In our case the symbol, defined by Eq. (38) above, is 
a step function jumping between A -|- 1 and A — 1, and 
the jumps occur at the zeros of A{9). We can assume 
that ^(0) > 0, as the local transformation Sj = b\ (which 
keeps the entanglement entropy invariant) yields A{0) — ^ 
—^(0). Using the notation for the zeros of A{9) by 9r, 
r = 1, 2, ... i? in an increasing order in the period (0, 27r] 
we can write the factors in (39) for the symbol (38) 



^(6) = (A + 1) 



A- 1 



<0 <2n + 9r, 



where 



^, = (-l)'^-Lln.^ + l 



27ri A - 1 

Indeed, one can easily check that the function given 
by (39) with the above defined ingredients has the value 
p{0) = A — 1 and alternates between A ± 1 with jumps 
at the zeros of A{9). Now, substituting our data in the 
statement of the conjecture we get the expression for the 
determinant 

det£>^(A) = (ln*)^L-«^' [] (l - e'^^^-^^^Y^ 

r=3 mod 2 

X n (l-e^^''-'^^y^\GB{l + P)GB{l-P)f . 



r^s mod 2 
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From this point, the calculation of the contour integral 
(37) is entirely identical to that of [18, 19], and the result 
for entropy asymptotics reads 

5i = |lnL-l J2 ln(l-e'(^-^'-))+ 

r=3 tnod 2 

i ln(l-e^(^-«^)) + |((l+7^;)-6/3ln2) , 

r^s mod 2 

(40) 

where the constants Te and I3 were given at the end of 
section III A. 



V. EXACT RESULTS FOR THE ENTROPY 
ASYMPTOTICS FOR CERTAIN 
NON-GAUGE-INVARIANT MODELS 

We now turn to discuss the cases of some non-gauge- 
invariant models. In the first two subsections we will 
determine the entropy asymptotics for chains that are 
Kramers- Wannier selfdual and for those that decouple 
to two independent majorana chains, by relating these 
cases to certain gauge- invariant models. In the last two 
subsections we will relate the entropy asymptotics of dif- 
ferent non-gauge-invariant models, by generalizing the 
XY-Ising transformation and doing local rotations. 

We will make use of the fact that one can write the 
general (translation-invariant) quasifrce Hamiltonian (4) 
in terms of the majorana operators defined by (7) in the 
following way: 

2N 

H = iJ2Tj,imjrrn (41) 

with the properties Tjj = —Tij e M and Tj+2n,i+2n = 
I for all n € Z. The transformation between the two 
descriptions reads 

T2j-i,2i-i = hm{Ajj + Bjj) 

T2j-i,2i = hie{^A,j+B,j) 
T2j,2i-i = hie{Ajj - Bj,i) . 

A. Kramers- Wannier selfdual models 

The Kramers- Wannier (or disorder) spin operators on 
a spin chain are defined in terms of the original spin op- 
erators (Pauli matrices) as 

I 



These spin operators also satisiy the Pauli commuta- 
tion relations = i'^^SkiSabc'^i- If the Hamil- 
tonian is invariant with respect to the above transfor- 
mation then it is said to be Kramers- Wannier selfdual. 
Such selfdual Hamiltonians always describe critical mod- 
els, an example is the critical point of the Ising model. 
A straightforward calculation shows that a quasifree 
Hamiltonian (41) is selfdual iff Tj^i = T^+i ;+i (recall, 
that translation invariance only implies Tj^i = Tj^2,i+2)- 
Or in other words, the selfdual models is the class of 
quasifree models, whose B matrix is real and the equality 
Kc{Ai,j + Bij) =Re(— Aij+i + Bij+i) is satisfied for ev- 
ery integers The two-point functions are then given 
by: 

(„,„,) = ,, + _L£e-'0-).;|||.« (43) 

with Te = J2n e"'''"To,n, thus in this case the block- 
Toeplitz matrix of the majorana expectation values re- 
duces to an ordinary Toeplitz matrix. Moreover, the 
above formula is (up to a factor of 2) identical to that of 
{b^jbi) for a gauge-invariant model with the symbol given 
as ^(6*) = iT{-e). 

Now, let us compare the calculations of the entropy 
from the matrix {niimj) in the general quasifree case (cf. 
Eq. (27)) and that from the matrix (blbj) in the gauge- 
invariant case Eq. (31). One can immediately conclude 
that the entanglement entropy of L spins in a Kramers- 
Wannier selfdual model defined by a matrix Tij equals 
half the entanglement entropy of 2L spins in a gauge- 
invariant quasifree model defined by Ai_j — iTi_j. Thus 
our result in Section IV applies to all Kramers- Wannier 
selfdual models as well. As an important example, we 
will apply this procedure in Section VI to obtain the an- 
alytic form of the entropy asymptotics for the critical 
Ising model with DM interaction. 



B. Directly decoupled majorana chains 

Next we turn to the case when the fermion chain 
decouples to two separate majorana chains. Prom the 
form (41) of the Hamiltonian we immediately sec that if 
T2i,2j+i = T2i+i,2j — for all which is equivalent 
to having purely imaginary matrices A and B, then the 
fermion chain decouples to two indc^pendent majorana 
chains: the one consisting of the odd modes and that of 
the even ones. Computing the symbol (24) corresponding 
to B"- = A" =0 gives 

/ -A<'(e)-B^(e) „ \ 

m = -i _^.(.HB-(.) , (44) 

V u Y^MW^wm ) 

hence the matrix is a direct sum of two Toeplitz ma- 
trices with symbols (— A° ± i3*)/| — A° ± | correspond- 
ing to the two uncoupled majorana chains. As in the 
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previous section, one can again relate the majorana ex- 
pectation values to the ground state expectation values 
of the bib^j operators of gauge-invariant models. Namely, 
we have 



{b]bl).A+B 
{b\bi).A-B 



-{m2jm2i) 

-(TO2j+im2(+l) 



(45) 
(46) 



where {■)-a±b stand for the expectation values in the 
gauge invariant models with H = J2jii~^ji ^ ^jiWjbi- 
Thus, by virtue of Equations (27) and (31) the entropy 
in the original model is given by the sum of entropies in 
the gauge invariant ones above. 



C. The generalized XY-Ising correspondence 

There is also a less direct way certain fermion chains 
can be decoupled into two independent chains. Suppose 
that matrix T in Eq. (41) satisfies the following proper- 
ties (for all 

Tii,ij-\ = 74i,4j_2 = T'4j_i^4j_3 = 0. 

By defining 



(1) -^(1) 



mo,_i =m. 



'2i-l 

(2) 
2i-l 



4i-3 ' "'■2i 



(1) - ^(1) 



m,/ = m 



: m. 



(2) 

Ai-2 > ^^21 



(2) 



m. 



4i ' 
(2) 

4i-l > 



(47) 



one can see that the original quasifree Hamiltonian with 
2N sites can be written as the sum of two other quasifree 
Hamiltonians with A'' sites: 



H 



2N 2N 

Here the components of matrices T^^' and T^^^ can be 
straightforwardly matched with the components of ma- 
trix T using the correspondence (47); it turns out that 
the decoupled subchains are also translation invariant: 

'^ij^'^'^^ = ^i+2j+2- This type of decoupling is the gener- 
alization of the famous XY-Ising correspondence [27] (for 
an other type of recent generalization of this correspon- 
dence, see [28]). 

Considering the ground state in the thermodynamic 
limit, this type of decoupling immediately implies that 
the entanglement entropy of 2L consecutive spins in the 
model defined by T equals the sum of the entropies of 
L spins for the models defined by T^^^ and T^'^K This 
method was used for deriving the entropy asymptotics of 
the critical Ising model from that of the critical XY chain 
[29], our result generalizes this. 



D. On general reflection- invariant models 

As we have discussed, even for refiection-invariant 

chains (Aij real, Bi j complex) there is no general for- 
mula for the entanglement entropy asymptotics. How- 
ever, as we mentioned in Section EIB, there is a formula 
for the saturation entropy in case the matrix B is real. In 
this subsection we show that a subclass of models with 
complex B can be transformed back to the real case. 

A transformation on the vector (m2j_i,m2j) = 
Uv^ with a constant matrix U G U{2) is called canoni- 
cal if the anticommutation relations {toj, to;} = 2Sji are 
preserved. For the two point functions it results in the ad- 
joint action {v^ ^v^) >-)■ U{v^ (^v^)U^ . Assume now that 
there are constants Cx,Cy,Cz G K with at least one of them 
non-vanishing, such that Ca;B"-{9)-CyA''{0)+C:,B''{9) = 
for all e [0, 27r). In this case there are rotations, 
which rotate the vector c = {c^, Cy, Cz) into c' = (0, 0, c^) 
(c'^ ^ 0) and consequently v = {B°- , —A'^ , B^ ) into 
v' = (B'",-A'^0) {as V ± c =^ v' ± d). Further- 
more, the Toeplitz symbol can be written as <^ = iM 1 — 
iP/V^Cl2^cF'^Va),{a = x,y,z), so the rotation can be 
done by the adjoint action of SU{2) on 2 x 2 traceless 
Hermitian matrices UrG{v)U^^ = G{Rv) {G{v) = G). 
This is exactly the above defined local transformation. 
Note, that the invariance of the entropy can be immedi- 
ately seen from the formula (36), which is invariant under 
the simultaneous transformation of all matrices by the 
adjoint action of any constant matrix (and the Wiener- 
Hopf factorization remains also valid). The general case, 
when A°',B°'^B'' are linearly independent Laurent poly- 
nomials of e'^, this method does not work. One could in 
principle try to follow a strategy similar to that of [21] 
as was done in [22] sketched in section IIIB. To obtain 
explicit results, where physical limits can be studied, is 
difiicult. 



VI. NEAREST NEIGHBOR COUPLING 

We will now look at the general quasifree model with 
nearest neighbour coupling and apply the above machin- 
ery to study its entanglement entropy. Our method yields 
analytic expression for the Ising model with DM interac- 
tion at the critical point, while for the general non-critical 
case we demonstrate that the ground state is not effected 
by the DM term, hence the results [22] apply. 

The Hamiltonian of the most general nearest neighbour 
spin chain that can be mapped to a quasi-free fermion 
chain is given by [36] 

H = ^((l + ^)ajaj+i + (l-7)a|a|+i + 



(48) 



The real parameters stand for the magnetic field h, 

the strength D of the DM current and the anisotropy 
7 G [0, 1]. The model is mapped by the Jordan- Wigner 
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transformation to the following fermionic one: 

J 

+l(b]b]^,-b,b,+,)-2hb]b,) . (49) 

One can analyse whether the one-particle spectrum de- 
termined from (17): 



A{9) 



^ Dsm6 + J (cos 6* - /i)2 + 72 gj^^ 



(50) 



vanishes or not at some to arrive at the following phase 
diagram: 



VII. SEEMING VIOLATIONS OF THE 
CALABRESE-CARDY FORMULAS 

In this section we will discuss two "anomalies" of the 
entropy asymptotics, which can appear at and in the 
vicinity of reflection-symmetry-breaking critical points 
and which seemingly do not fit the Calabrese-Cardy for- 
mulas. The first concerns the growth of the saturation 
entropy as we approach such critical points, while the sec- 
ond is about the breaking of refiection symmetry in the 
finite-size scaling of the entanglement entropy. We will 
discuss how we can interpret these anomalies to keep the 
validity of the Calabrese-Cardy formulas. 




The parameter D' is defined by 



D' = y/D^ + 1-72 



(51) 



and the critical regions are (i) the connected one between 
the D'^ = h? paraboles and the D'"^ = 1 hue and (ii) 
the \h\ — 1 line segments. One immediately observes 
from the form of the two point correlations (24) that the 
ground state of the non-critical regions are given by the 
XY model as M(6') = for all 6^ € [0, 2it) and it is via 
M{9) that {rrijmi) depends on D. 

Let us turn now to the special case of 7 = 1 that is, 
the Ising model with DM term. The phase diagram is the 
same as above with D'^ = D^. The case h — I can be 
solved by noticing that this case belongs to the class of 
models with T matrix of the genuine Toeplitz type. We 
have only two non-zero elements Ti = 1/2, T2 = — Z)/4, 
which gives iT{9) = 4 sin 0(1 — Dcos9) for the numer- 
ator of the symbol (43). The entropy asymptotics for 
this class is given via Sl ~ S'^*"^/2, as explained af- 
ter Eq. (43), the superscript refers to the gauge-invariant 
model, whose Toeplitz symbol is A. Its entropy is given 
given by Eq. (40). The final result reads 



Sl - 
with 

CisDM 



ilnL+^ln(l-^ 



ilnL + C/, 



CisDM \D\ > 1 
Il?l< 1 



(52) 



- (1 +7B + (1-6/3) ln2) 



0.726067 



Cis = -(l+7£ + (2-6/3)ln2) « 0.478558 
6 



A. Anomalous behavior of the saturation entropy 

As mentioned in the Introduction, the formula for the 
saturation value of the entanglement entropy of a block 
of spins near a critical point reads 



Ssat = 3 



const , 



(53) 



where c is the central charge belonging to the critical 
point. We have seen in the previous section that con- 
sidering the region < I?' < 1 (see Eq. (51) for the 
definition of D'), the XY model with DM interaction 
is critical when h = ±1 and the corresponding central 
charge is cxy-dm = 1, while for the model without DM 
interaction {D = 0) the central charge of the critical line 
(at h — ±1) is cxY = 5- However, we have also shown 
that when /i ^ ±1 the ground state does not the depend 
on D in the non-critical region < I?' < 1. Approaching 
the critical h — ±1 line in this region, the divergence of 
the saturation entropy (which is hence independent of D) 
is consistent with formula (53) in case the central charge 
is c = i, as can be seen from the results in [22, 30]. Hence 
the formula is not valid for the XY model with DM in- 
teraction, since for that model the central charge is 1. 
This situation is typical for quadratic models with reflec- 
tion symmetry breaking: We have seen in Section HA 
that the ground state of the Hamiltonian Eq. (4) and the 
central charge do not depend on ImAij at a non-critical 
point, while they may depend on it at a critical one. 

Hence, in order to understand the failure of formula 
(53), and to formulate a possible reinterpretation in the 
case at hand, let us first look at an other anomaly, which 
is, in some sense, similar. In the XX-modcl with trans- 
verse magnetic field 

i 

there is quantum phase transition at the points h = ±1. 
In the region —l<h<l the ground state of the 
model is critical with algebraically decaying truncated 
correlation function C^^{n) = (crfaf^_„) — Cq^ (where 

= hm, 
totics Sl 



erf erf ^n)) ^ diverging entropy asymp- 
ilnL -|- k, while outside this region the 
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ground state is cither the aU-spin-up or all-spin-down 
state (depending on the sign of h). Hence approaching 
the critical region from the non-critical one, the satu- 
ration entropy will not diverge, however, this does not 
contradict formula (53), since there is no diverging cor- 
relation length either. When we enter the critical re- 
gion the state will change suddenly in such a way that 
the correlation function C^^{n) that was identically zero 
in the noncritical region suddenly will be nonzero and 
even quasi-long-ranged (decay algebraically with n), i.e., 
an "infinite correlation length" appears instantaneously. 
This is, in some sense, a degenerate situation, because 
considering a bigger parameter space, e.g., the XY model 
with transverse magnetic field 

H = Y.{.1+ 7)«+i + (1 - 7)^M'+i + hat , 

i 

and approaching the critical line (— l</i<l,7 = 0)by 
fixing the value of h (between 1 and —1) and taking 7 
1, we will observe a diverging correlation length and a 
diverging saturation entropy satisfying formula Eq. (53), 
as can be seen from the results in [22, 30]. 

A similar situation, h\it in a more complicated form, 
occiirs in the XY model with DM interactions. At 
the critical point, the correlation functions C^^(n) and 
C''(n) = {crfo-f^n ^ ^i^i+n) decay algebraically. How- 
ever, away from criticality C"^ (n) is identically zero, while 
C^^{n) behaves "in a normal way", i.e., it decays ex- 
ponentially with n and the correlation length length, 
characterizing the exponential decay, diverges as one ap- 
proaches the critical line h ^ 1. We can think that there 
are two independent critical "modes" both with c = ^, 
one is behaving normally, the other in an anomalous way, 
hence in Eq. (53) we should only insert the central charge 
of the normally behaving mode. This picture could be 
made more convincing and precise, if one could show 
that, similarly to the previously mentioned XX case, this 
anomalous behaviour is a degenerate one by considering a 
bigger subspace, e.g., the XYZ chain with magnetic field 
and DM interaction, i.e., adding the term ^ . Aaf cr?^^ to 
the Eq. (48). We conjecture that approaching the men- 
tioned critical point in such a bigger parameter space the 
generic behavior induces normally diverging correlation 
lengths for both C^^(n) and C'^(n), denoted by and 
^j, respectively; and the entropy will scale according to 
a "generalized" Calabrese-Cardy formula of the form of 

S'l = ^ In ^xx + ^ In G + const . (54) 

We have started to study this conjecture numerically, and 
the results will be the subject of a forthcoming publica- 
tion. 

B. Breaking of reflection invariance in the finite 
size scaling of the entanglement entropy 

The other feature we will investigate is whether (and 
to which extent) the breaking of reflection invariance can 



be observed as a finite size effect in the scaling of the 
entropy. More precisely, consider a finite spin chain with 
a quadratic Hamiltonian of Icingth N, and compute the 
entropy S{L, N) of the restriction of the ground state to 
the first L sites. The Calabrese-Cardy formula (3), which 
has been confirmed (up to subleading corrections) ana- 
lytically and numerically for many reflection-symmetric 
models [6, 7, 31], suggests a reflection-invariant form 
S{L,N) = S{N - L,N). The question we ask is 
whether this symmetry of the ground state can be bro- 
ken for a quadratic Hamiltonian, which is not invariant 
and whether the symmetry breaking survives the limit 
N ^00 (with L/N fixed)? 

First, we should notice that the reflection invariance of 
the entropy function S{L, N) can only be broken if nei- 
ther the matrix A nor B is real for the following reasons. 
We saw that the Hamilton operator (49) is invariant (and 
so is the unique ground state) unless ImA ^ 0. For the 
case ImB = 0, one should consider the transformations 
bi — > b]\r-i and hi — >■ &J and determine the; transformed 
density matrices restricted to the first L sites of the chain. 
For the case ImB = they are identical (both transfor- 
mations lead to changes A ^ A''' , B ^ — i? in the Hamil- 
tonian) . The first one corresponds to the reflection we are 
interested in, whereas the second is a local transforma- 
tion of the chain and those preserve the entanglement 
entropy. 

As noted in Sec. VI, the nearest neighbor quasifree 
Hamiltonians can always be transformed by local trans- 
formations such that B is real. Hence to have a symmetry 
breaking entropy function we have to consider next-to- 
nearest neighbour Hamiltonians. The particular Hamil- 
tonian we investigated was 

JV 

i=l 

with the following parameters: ti ~ 7 + 28i, t2 = 4: + 5i, 
pi = ll + lOi, p2 = 3 + 4i, and h = 12. 

The numerical results depicted in Fig. 2 demonstrate 
that the reflection symmetry of the entropy function 
S{L, N) is indeed broken. However, it is also visible that 
the deviation 5(i, N) — S{N — L, N) goes to zero in the 
limit — >■ 00 for any fixed L. Moreover, we can see 
in Fig. 3 that in this limit the curves nicely converge to 
the Calabrese-Cardy formula. Hence we conclude that, 
according to our numerical results, the reflection symme- 
try of the entropy function can be broken, but its scaling 
limit shows no such breaking, and the Calabrese-Cardy 
formula is valid. 



VIII. SUMMARY AND OVERVIEW 

In this paper we studied the entanglement entropy 
asymptotics of spin chains that can be mapped to 
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FIG. 2: The figures show AS' = S{L,N) - S{N - L,N), 
the difi'erence between the entropies of the ground state re- 
stricted to the left- and rightmost L modes of the chain 
with N = 1024, 2048, 4096. The black lines indicate AS = 
±10~^. The plots show, in all three cases, the range x = 
L/N G [0.04,0.96]. Beside the feature that the function 
A'' I— >■ AS(A'') decreases for any fixed L, it also has an ad- 
ditional oscillating structure. One finds that the analytic 
function p{N/L + N/{N - L)) cos 27rL(A/Ar - 1/5) fits this 
structure rather well with a suitable constant p for the ampli- 
tude and the A'^— independent wavelength A. The explanation 
of this behaviour is under investigation. 




0.1 0.2 03 OA 0.5 0.6 0.7 0.8 0.9 



FIG. 3: Here we can see the entropies S{L, N) for TV = 
256, 512, 1024, 2048, 4096 and the corresponding Calabrese- 
Cardy curves. The central charge is a fit parameter and it 
converges to the physical value c = 1 fast (the deviation de- 
creases roughly linearly with N from 0.038 at iV = 256 to 
0.003 at TV = 4096). 

quasifree fermionic models given by the sum of a gauge- 
invariant term (parametrized by a selfadjoint matrix A) 
and a non-gauge-invariant one (parametrized by an anti- 
symmetric matrix B). Many models of physical impor- 
tance belong to the class of complex A (and B), implying 



the breaking of reflection symmetry). The entanglement 
properties of these systems have hardly been addressed 
in the literature before, hence we concentrated on these 
cases. 

We have determined the two-point functions of the ma- 
jorana operators in complete generality. A novelty follow- 
ing from this investigation is that the ground state can 
only be reflection symmetry breaking if it is critical. 

We have been able to write down the analytic ex- 
pression of the entropy asymptotics for the most gen- 
eral gauge-invariant models, and also extended these re- 
sults for certain non-gauge-invariant models. A detailed 
investigation of the nearest neighbor case was carried 
out. We have derived the explicit form of the entan- 
glement entropy asymptotics for the Ising model with 
Dzyaloshinskii-Moriya interaction at the critical point, 
which was unknown until now. In the noncritical regime, 
we demonstrated that the ground state is independent of 
the DM coupling, thus the entropy asymptotics given in 
[22] without the DM term is valid also here. This indi- 
cated violations of the formula for the saturation entropy 
^ clogf near the critical point \h\ — 1. We have given a 
possible "physical explanation" for this. 

Concerning the general landscape of the block entropy 
asymptotics of quasifree models, we extended the gen- 
eral knowledge to a large extent, nevertheless, the gen- 
eral case remains to be a difficult unsolved mathematical 
problem [37] . Even when specifying the discussion to the 
nearest neighbor case, there remains a surprisingly large 
region of the critical regime, for which the scaling of the 
block entropy still remains an open problem. 

Finally, we carried out numerical checks for the inves- 
tigation of finite size effects. We used a model Hamil- 
tonian with next-to-nearest neighbor interaction, which 
exhibited reflection symmetry breaking in the finite-size- 
scaling of the entanglement entropy. The deviation was 
demonstrated to converge to zero quickly by increasing 
the size of the chain, while the block entropy converged 
to the asymptotic Calabrese-Cardy formula. 
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